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Abstract 



We investigate the issue of charge neutrality of the CFLK phase of dense 
quark matter using the low energy effective theory of high density QCD. We 
show that the local electric and color charge neutrality of the ground state 
in a homogeneous color superconducting medium follows from its dynamics. 
We also consider the situation of a spatially inhomogeneous medium, such 
as may be found in a neutron star core. We find that spatial inhomogeneity 
results in the generation of electric fields, and positrons and/or electrons may 
be present in the ground state. We estimate the concentration of charged 
leptons in the ground state to be n e ~ 10 2 cm~ 3 and consider their influence 
on the opacity of the medium with respect to the modified photons. 

I. INTRODUCTION 

It is now understood that at asymptotically large baryon number densities, the ground 
state of 3- flavor, 3-color massless QCD is the Color-Flavor Locked (CFL) phase [1,2]. (See 
[3] for comprehensive reviews.) In this phase, all 9 species of quarks with momenta close to 
the Fermi surface undergo BCS-like pairing. The resulting quark-quark condensate spon- 
taneously breaks the original SU(3)c X SU(3)l x SU(3)r symmetry down to the diagonal 
subgroup, SU(3)c+l+r, causing the gauge bosons of the original SU(3)c group to become 
massive, and also breaks the U(1)b symmetry associated with conserved baryon number. 1 
One linear combination of gluon and photon remains massless; the associated unbroken 
gauge symmetry is referred to as U(1)q, denoting the "modified" electromagnetism. 

A nonzero mass for the strange quark encourages the system to reduce the strange quark 
number density relative to the density of up- and down-type quarks. It has been argued 
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X A local gauge symmetry cannot really be spontaneously broken [4]. We have to fix a gauge to 
define expectation values of gauge variant quantities. 
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that the system responds to this stress by forming a kaon condensate in the ground state 
[5-7]. Kaon condensation allows the strange quark number density to be decreased without 
the costly breaking of pairs in the CFL background. 

In particular, it was found that in the presence of small enough chemical potentials for 
the electric and lepton charges K° condensate was present in the ground state (the CFLfT 
phase) [6,7]. The CFL and CFLK phases are distinct, as hypercharge is spontaneously 
broken in the latter, but not in the former. Some recent reviews are listed in [8]. 

Recently, some attention has been devoted to the issue of realization of charge neutrality 
in the high baryon density matter. First, it was argued that for a macroscopic system the 
difference between being in a color singlet state and being in a state with equal number den- 
sities of red, blue and green color charges is negligible [9]. Alford and Rajagopal constructed 
an expression for the free energy that incorporated general features of the CFL phase [10], 
while Steiner, Reddy and Prakash used NJL model supplemented by diquark and the t'Hooft 
six-fermion interactions [11]. Neither approach included gauge bosons as dynamical fields. 
Instead, to ensure charge neutrality of the CFL phase, chemical potentials // 3 and /i$ were 
introduced for the charges corresponding to T3 and T§ generators of the SU(3) C group, along 
with hq, a chemical potential for the electric charge. In both papers it was found that /X3, 
/is and fiQ had to be adjusted to ensure local charge neutrality of the bulk quark matter. 

Motivated by the ideas of [10], in this paper we employ the low energy effective theory 
of QCD at high baryon number density developed in [12,13,6] to address the issue of charge 
neutrality in the high baryon density matter. In particular, in this approach gauge fields 
are treated as dynamical degrees of freedom. One expects that in a color superconducting 
state, color charge neutrality will follow from the low energy dynamics of the system, and 
that is what we find for a homogeneous case. 

In section 2 we show that this is indeed the case for a homogeneous medium; in section 3 
we present the analysis for the case of an inhomogeneous medium, such as one would expect 
in a compact star, and we find a surprising result that nonzero electric fields exist in the 
bulk in this case. 

II. CHARGE NEUTRALITY IN HOMOGENEOUS MEDIUM 

QCD at asymptotically high baryon density has several effective field theories valid for 
different ranges of excitation energies. For energies well below the quark number chemical 
potential \x heavy anti-quarks (mass ~ 2/x) may be integrated out in a systematic manner 
leading to the High Density Effective Theory (HDET) developed in [15,16,6]. 

For energies smaller than the superconducting gap A quasi-particle and quasi-hole states 
in the vicinity of the Fermi surface may, in turn, be integrated out leaving the Nambu- 
Goldstone bosons due to breaking of various flavor symmetries in the ground state as the 
relevant degrees of freedom [12,13,6,17]. For instance, in the CFL phase we have an octet 
due to breaking of the chiral symmetry and a singlet due to U(1)b breaking. See [18] for 
a review of the subject. The low energy effective theory of dense QCD is similar in spirit to 
the effective theory of an ordinary (QED) superconductor (Ginzburg-Landau theory). See, 
for example, a pedagogical introduction by S. Weinberg [14]. 

The leading terms of the low energy effective Lagrangian of QCD at high baryon number 
density are 
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C = -G- Tr ([X^D^X) 2 - h^X^DfXf^j + Tr {{Y ] D*Y) 2 - h^DfY) 

-\F, U F^ -±Tr (G„ u Gn, (1) 

where X and Y are the 3x3 unitary matrix valued composite fields describing oscillations 
of the left- and right handed quark-quark condensates about the CFL ground state. Under 
the original symmetry group SU (3) c x SU (3) L x SU(3) R , X transforms as (3, 3, 1) and Y as 
(3, 1,3). The pion decay constant f w ~ 0.209^ was computed in [13]. The gauge covariant 
derivatives are 

D L V X = d u X + ig s X{A c v ) T + ieA e v m QX + i8 u0 ^-X, 

D*Y = d v Y + ig s Y(At) T + teA e u m QY + t5 u0 ^-Y, (2) 

where A c v and Al m are the gluon and photon gauge fields, respectively; A c v = A c au t a , where t a 
is an SU(3) C generator in the fundamental representation. We treat quark mass matrix M 
as an external (spurion) field that transforms as (1,3,3). It was arg ued in [6] that M 2 /2fx 
terms should appear in the covariant derivatives (2), and, thus, that in the low energy 
effective theory m 2 s /2fi plays the role of chemical potential for strangeness. These are the 
terms responsible for the formation of kaon condensate in the ground state. In this paper 
we neglect up and down quark masses and set 

M = diag(0,0,m s ). (3) 

The quark electric charge matrix is Q = diag(2/3, —1/3, —1/3). 

A generic term in the effective Lagrangian is made of various SU (3)c x SU (3) l x SU(3)r x 
U(1) b invariant combinations of powers of f 2 A 2 Y^(^) k Y, f 2 A 2 X^(^) n X, G„ v and F^. 
The terms involving both X and Y and the ones containing the field strengths are further 
suppressed by powers of g s and/or e [13]. 

It is more convenient to rewrite the Lagrangian in the following form [12] 



C = C XYA + £ s , (4) 



where 



Cxya = - yTr (x*X + yty + ^(^0- + eA™Q)X + iY\^~ + eA^Q)Y + 2ig a {A^ 
+ ^I Tr (x^VtX + Y^ViY + ieAt m (X^QX + Y^QY) + 2ig s (A c i ) T ) 2 
-\F„ U F^ -±Tr (G„ u Gn, (5) 

and 

£ s = Jl Tl (W - Y^Y + - iY ^-Y + ieA$(X*QX - Y^QY)) " 

4 V 2^(x 2^tx / 
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Here 



and 



+ ^Tr (xWiX - Y ] ViY + ieA®(X ] QX - Y*QY) 
f 2 

^Tr (s + iWJl - iEW R ) (£ f + iW R T) - itfW L ) 

-S-^Tr (V.E + ieAf[Q, £]) (v^ + i&4?[Q, £ f ]) . (6) 



^ = e~4Q + ^ (7) 



^ = e ~4g + ^, (8) 



2/i ' 

is the Uq(1) gauge field with the coupling constant 

e = fit . (9) 
V^ 2 + 4e2 

We have introduced the color singlet field S = XY^ = exp[2in a t a / f n ], where 7r a 's are 
the pseudo-scalar octet of Nambu-Goldstone bosons which arise from the breaking of chiral 
symmetry. E transforms as (1,3,3). We neglect the Ub(X) Nambu-Goldstone boson terms 
in the Lagrangian as they are irrelevant for the gauge charge neutrality. 

We see that Cs only depends on light fields, while Cxya has both light and heavy degrees 
of freedom in it. The convenience of the representation (4) is due to the fact that at low 
energy the heavy gauge fields may be integrated out. Then Cxya piece vanishes leaving us 
with £s only. Now one may proceed to solve equations of motion to determine values of 
various (classical) fields in the ground state. It was found in [6,7] that in the ground state 
E field takes on a non trivial value 



K o = i \ (10) 




corresponding to CFLfT phase. We note that in the approximation m u = rrid = the K° 
condensate in the ground state is maximal and formed for arbitrarily small value of m s [7]. 

From Cxya we see that Xi^jf-X + Y^^fY acts as a source term for the temporal 
components of gauge fields. This follows from M 2 /2/i terms being part of the covariant 
derivatives (2). We may also see this from the matching calculation of the coefficient of the 
operator 

9s A c a Tr ((t a ) T M) , (11) 

where 
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K vD_0_00_Q 



FIG. 1. This is the HDET diagram that produces operator (11) at one loop level. Fermion lines 
with triangle insertions denote the anomalous quasi particle propagators. Hollow circle stands for 
insertion of M defined in (12). Normal one loop diagram contribution vanishes. 



M = x^ x + r f y, (12) 

in the low energy effective Lagrangian (4). The calculation yields a non zero coefficient of 
/ 2 /4 that comes from the anomalous contribution to the one-loop tadpole diagram with one 
insertion of A4 as shown in Fig. 1. The diagram is superficially proportional to A 2 , but the 
loop integration is infrared sensitive and picks up a compensating 1/A 2 . 

We note that spatial components of gauge fields do not have sources and vanish in the 
ground state. 

To proceed it is convenient to define a unitary matrix £, such that £ 2 = S, and an SU(3) 
matrix U . Then X and Y may be presented as 



X = £U, 

Y = ?U (13) 



We fix SU(3)c gauge by setting 



U=l, (14) 
which corresponds to the unitary gauge [12]. Then we have 



( V2 

1 i | (15) 
V % 1 



and X = £ K o, Y = £^ . 

In the homogeneous medium all fields are constants and the Lagrangian reduces to 




U° + ^°^7-^ ot + 2 9s(A c f + 2eA™Qj (16) 



Now let us proceed to state the main result of the paper. The charge densities are defined 

as 



J " = ~w c - (17) 

while the equations of motion for the temporal components of the gauge fields in the homo- 
geneous case read 

-C EE —r-C = 0. (18) 



SA dA 

Then color and electromagnetic charge densities vanish in the ground state since the time 
components of the (classical) gauge fields must satisfy their equations of motion. Note that 
vanishing of charge densities in the ground state is a gauge independent result and does not 
depend on the particular (gauge dependent) values of the gauge fields in the ground state. 

From the Lagrangian (16) we see that the charge density due to the nonzero strange 
quark mass (the first two terms in (16)) induces non zero A fields in the ground state. 
So, when substituted into the expression for the charge density, the Debye screening term 
(the last two terms in (16)) and the m s -induced charge density term cancel each other 
and the corresponding charge density J vanishes. This is just what we expect to see in 
a superconductor. Qualitatively, this is similar to the picture from [10]. There the linear 
in the chemical potentials charge density terms (analogous to the Debye screening terms in 
(16)) cancel the m s -induced quark charge density pieces in the expressions for the charge 
densities thus ensuring charge neutrality. 

Now let us solve equations of motion for the A 's. The formulas for the gauge fields will 
be needed in the next section. Let us note first that lagrangian (5) (and therefore (16)) does 

not depend on the A®. This is a manifestation of the fact that the medium is Ug(l) neutral. 
The equation of motion for A® is V 2 A® = and the solution is 

4 = C, (19) 

where C is a constant. Then solving the equations of motion $j^jC = 0, with C from (16), 
yields for A 's (A = 0) in the original basis 



em s Ce 
3 4/i^e 2 g s ' 

e 2 m 2 Ce 



03 = 
08 = 



4v / 3/ifi' s e 2 V3g s ' 
e 2 m s 2 Ce 



S/j,g s 2 e e 



+ — . (20) 



All other gauge fields vanish. In the pure quark matter Aq may take an arbitrary value. 
However, when one takes into account electrons, the requirement of CFLK electric charge 
neutrality gives 2 



2 I am grateful to M.Alford and K.Rajagopal for pointing out this fact to me. 
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so that <p em = 0. The remaining fields take on the following values 



^ = _^L_. (22) 

Finally, let us make an important observation that due to the constraint imposed by 
the gauge invariance on the form of operators that are allowed to appear in the low energy 
effective theory Lagrangian, the expressions in (20) are not corrected by the neglected higher 
order terms in the effective Lagrangian (1). 



III. SPATIAL INHOMOGENEITY AS A MECHANISM FOR GENERATING 
ELECTRIC FIELDS IN HIGH DENSITY QCD 

Our formalism also allows us to consider the case of inhomogeneous medium such as may 
be found in a neutron star core. We assume \i to be an unspecified, spherically symmetric 
function /z(r) with some typical size of inhomogeneity L. We use L = 1 km for the numerical 
estimates, which is the typical length scale of a neutron star core. The characteristic length 
scale of the ground state is the screening Debye length \ D ~ {gsl^) 1 ~ Ifm. Since Xd/L ~ 
1CT 18 , the system is only very weakly inhomogeneous. We will analyze the system in an 
expansion in e = Xd/L. For now we will neglect the electrons and consider only the quark 
matter. 

Let us state the main result right away. To leading order in e, the equations of motion 
for the gauge fields are = 0, with C given by (16), but with /j, replaced by /i(r). The 

solutions are given by (20) with ji replaced by /i(r). Then 

(f) ~ m 2 s fix (r ) + const + £>(e 2 ). (23) 

We immediately notice that, for example, 

o 

ID 

E r = -d r <f) em MO. (24) 

jiL 

So, we make a surprising observation that the spatial inhomogeneity gives rise to non zero 
electric fields in the bulk of quark matter. 

It is useful to consider formation of such weakly inhomogeneous state, that is a situation 
where \i = /i(t, r) evolves as schematically shown in Fig. 2. Starting from a homogeneous 
density profile at t — 0, the system grows denser in the central region as time elapses. We 
set the typical time evolution scale T to be much bigger than any relevant QCD time scale, 
the typical size of the inhomogeneity L is much greater than as discussed above. 

Let us consider, for example, the equations of motion for A e ™. We have 




r 



FIG. 2. Starting from a homogeneous profile at t\ = density increases in the central region 
while remaining constant far away from the center. By the time t±^>T density profile approaches 
its stationary asymptotic form /i(r) and j r defined in (30) vanishes, but the charge density, (28), 
remains non zero. 



dF v» = C, (25) 

f 1 A A p.m. ' v ' 



Q A em 

where C is given by Eq. (4). We seek solutions in the form 

em = 0oe° + 2 e 2 + ..., 

A em = A e° + A 2 e 2 + .... (26) 

To leading order in e the equations of motion for the gauge fields are -^-C = 0, with C given 
by (16), but with \x replaced by fj,(t,r). The temporal components of the fields are given by 
(20), and we find 

m 2 s e 2 ^ Ce 
° 3egifi(t,r) e ' 
A) = (27) 

where constant C is denned in (19). At order e 2 we find the following expressions for the 
electromagnetic current components 

J° em = -V 2 0o (28) 

and 

j em = cW0 o . (29) 

Note that d^j£ m = 0, so that the electric charge conservation is satisfied. Thus, there exists 
a current with nonzero radial component 

2-2 2 
777 P 777 

*<'•'> = ~ sir (30) 
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Now we see that formation of an inhomogeneous density state is accompanied by the po- 
larization of the system; the positive charges move outward and the negative ones - inward. 
Let us emphasize that during this process the system remains charge neutral. As the density 
profile reaches its stationary form, j r (t, r) vanishes but the charge density, (28), remains non 
zero. 

Let us consider a stationary inhomogeneous case \i = /z(r) such as the asymptotic density 
profile fi(r, £4) in Fig. 2. So, we consider a spherically symmetric inhomogeneous region of 
size R> L embedded into infinite homogeneous CFLf^ medium. 

As we have mentioned earlier, a spatial inhomogeneity of the density generates electric 
fields in the quark matter. Then it is possible that the light charged leptons, i.e. positrons 
and electrons, may be attracted into the system to neutralize the U em (l) electric field created 
by the quark matter. To address this question we add 

C e = 4>{iYd v - m e )^ - e^YMT (31) 

to the Lagrangian (4), where vp is the usual electron-positron Dirac field and m e is the 
electron/positron mass. Now, using equations of motion we integrate out all the remaining 
heavy gauge fields. We are left with the low energy effective Lagrangian 

C = ^i 1 v d u -m e )^- l -A%V 2 A<l 

( *o~ m 2 e 2 \ 2e 4 n 2 
-n- e A$e + - -— 32 

where n s = ip^ip is the fermion density. We notice that fermions appear in the Lagrangian 
(32) as coupled to a gauge field with nonzero temporal component equal to 

Mr) = Al(T) + SFY (33) 

One may think of the fermions as being in the presence of an effective chemical potential 
— eA. Then we may assume that in the ground state fermions form a Fermi sphere with 
Fermi momentum kp being a slow function of r (Thomas- Fermi approximation). In the 
uniform case the electric charge neutrality of CFLK requires that no electrons should be 
present in the ground state meaning that 

A = A$ + ^ = 0, (34) 

which is consistent with (21). 3 

We find the Fermi momentum from the equality of energy of a particle on the top of the 
Fermi sphere inside the system at a given point and a particle at rest far away from center 
at some r = R, where the density is constant (uniform density region far away from the 
center in Fig. 2). The equation is 



3 Strictly speaking, the condition is \eA\ < m e . We neglect m e compared to m 2 //i ~ 50MeV. 
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y 'k F {r) 2 + m 2 e + eA{r) =m e + eA(R) = m e . (35) 
Then in Thomas- Fermi approximation, ignoring 0(n 2 ) term in (32), we have 

n e (r) = 7^k F (rf = ^((eA(r) - m e f - m 2 e f . (36) 
Solving the equation of motion for A® 

V 2 4 = -en e (r) = ~((4(r)e + - m e ) 2 - m 2 e f' 2 (37) 

by the expansion in e yields 

A$ 
n e 

Since V 2 (l/yu(r)) ~ l/L 2 /i(r), we have n e ~ m 2 s /L 2 n ~ 10 2 cm~ 3 and k Fe + ~ 10~ 7 MeV 
for ji = 400MeV and m s = 150MeV. An a posteriori estimate gives ~ m e m 2 /L 2 fi for the 
size of the neglected n\ term, which is subleading in \d/L compared to the size of the 
n s eA(r) term in (32) which we retained. 

From (38) we can see that for the case of our toy model density profile the total electric 
charge enclosed in the inhomogeneous region vanishes by Gauss' law and we don't have to 
worry about the energy contribution from the long range electric fields. 

We estimate mean free path I of Q photons propagating in such a medium. For the case of 
low energy incoherent scattering of photons off of e + we get / = l/n e a ~ m 2 /a em n e ~ 10 19 cm 
which is much larger than the radius of a star. So, the medium is still transparent. Let us 
stress, however, that taking into account the spatial inhomogeneity results in a qualitative 
change in the ground state. Namely, light Q charged particles (positrons and/or electrons) 
are present in the ground state, while the homogeneous CFLK phase is a perfect Q insulator. 



m 2 e „„. m7 3 



+ ^( T77CT- ) 



3g 2 fj,(r) L 4 / 3 yU 2 / 3 m 
^ 2 %^ + ° ( LWV/3 me )- ( 38 ) 
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Postscript 

After this work had been completed and the paper was being prepared for submission I 
became aware of Ref. [19], where a similar conclusion about the charge neutrality of the 
color superconducting quark matter was reached. 
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